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ABSTRACT 


Direct application of two-phase heat transfer in the liquid cooling of clectronic 
components in flouorinated hydrocarbons (FC-72), is severely inhibited by the excessive 
amount of superheat required to initiate nucleate boiling. To conduct an experimental 
study of nucleate pool boiling of F'C-72, an experimental test chamber was constructcd. 
This chamber utilized five horizontal platinum wires of 0.05 mm diameter spaced 2.0 
cm vertically from cach other. The lowest wire was progressively heated from the na- 
tural convection region through nucleate boiling, and a study was made on the effects 


of the boiling wake plume on the heat transfer rate of the upper wires. 
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NOMENCLATURE 


Description 

Surface Area of Test Surface 
Thermal Expansion Coefficient 
Fluid’s Specific Heat 

Wire Diameter 

Acceleration due to gravity 
Heat Transfer Coefficient 
Thermal Conductivity 

Wire Lengths 

Heat Transfer 

Heat Flux 

Resistance 

Temperature 

Rayleigh Number 

Prandlt Number 

Nusselt Number 


SUBSCRIPTS 


liquid 

vapor 

Saturation Conditions 
Subcooled Conditions 
Bulk Fluid Conditions 
Platinum Wire 

22 Resistor 

Surface 


GREEK SYMBOLS 


Dynamic Viscosity 

Kinematic Viscosity 

Surface Tension 

Fluid Density 

Uncertainity in measurements. 


I. INTRODUCTION 


A. HISTORICAL 

Ever since the development of the first practical solid state integrated circuit in 1959, 
there has been a continuous emphasis on packaging an ever higher density of logic cir- 
cuits on a single substrate. Corresponding with this increase packaging density is an 
increasing power dissipation requirement along with an increase difficulty in the thermal 
management capabilities. Initially, thermal management techniques consisted of noth- 
ing more than a heat sink and natural convection air cooling. This naturally progressed 
into heat sinks with forced convection air cooling, and with the construction of large 
mainframe computers and first generation super-computers the use of indirect liquid 
cooling was initiated. The next practical method for heat removal is direct liquid cooling 
or immersion cooling of microelectronic components using dielectric liquids. 

Although direct liquid cooling has been successfully employed for nearly 40 years in 
large hi-powered microwave equipment, klystron tubes, travelling-wave guide tubes, and 
power transistors [Ref. 1 : p.12], only recently has it been employed for use in the ther- 
mal management of microelectronics. The first computer systems to successfully employ 
immersion cooling of microelectronics were the CRAY-2 and the ETA-10 supercom- 
puter [Ref: 2: p. 8]. Today the research in direct immersion cooling is continuing with 
particular emphasis on utilizing dielectric liquids as a cooling medium, in particular the 
prefluorinated FC liquid family of the 3M corporation [Ref. 3]. 


B. BACKGROUND 

The advantages offered by direct immersion cooling is self evident from Figure 1 
By utilizing forced convection direct liquid cooling with FC-77, the CRAY-2 has been 
able to reduce it’s required volume by an order of magnitude as compared to it’s indirect 
liquid cooled competitors. This was being achieved in the CRAY-2 while at the same 


time maintaining the same power dissipation of the competition, [Ref. 1]. 


Mainframe, 
Supercomputer 


CRAY-2 @ 
(immersion cooled) 


Electric power consumption 


Office Computer 





0.1 10 100 


Volume m2 
Figure 1. Power consumption vs Volume: Of present computers 


As previously stated, the CRAY-2’s packaging density was achievable due to direct 
liquid immersion cooling, however, it did use a forced convection cooling scheme. A 
thermal management scheme that relies on boiling heat transfer would provide an even 
more efficient heat transfer capabilitv, and has been the attention of considerable re- 
search in recent years. Unfortunately, these prefluorinated dielectric fluids possesses 
anomalies which have held back their widespread use and acceptance of direct immersion 
cooling. 

In comparison to water and many other coolants, these prefluorinated dielectric 
liquids posses a lower value for their thermal conductivity, and heat of vaporization, 
they also have an extremely small surface tension which translates into a near Zero wet- 
ting angle and therefore, a high wettability on most surfaces. This last feature results in 
the elimination of many potential nucleation sites which in turn results in a incipience 
super heat excursion which results in a very high wall temperature. 

In this study, heat transfer enhancement, and therefore a surface temperature re- 


duction, due to bubble pumping in FC-72 will be investigated. 


C. PREVIOUS WORK 
An extensive amount of research has been conducted in the boiling regime of 


dielectric liquids. The most helpful of the literature used will be noted. 


You, Simon and Bar-Cohen [Ref. 4 and Ref. 5 ] conducted an analyses on the ef- 
fects of pressure, subcooling, and dissolved gas content on the boiling heat transfer, and 
the boiling incipience of dielectric fluids. In their experiment they used FC-72 as the 
working fluid and a 0.lum thin film of platinum on a 0.51mm diameter quartz cylinder, 
with the platinum film acting as both heater and temperature sensor. They concluded 
that: 


1. An increase in pressure of 56kPa reduces wall superheat by approximately 3°C in 
the fully developed nucleate boiling region. 


2. Subcooling has little effect on the fully developed nucleate boiling region. 


3. Nucleate boiling hysteresis is observed for both the subcooled and gas-saturated 
cases. The hysteresis decreases as the gas content increases. 


4. Incipience bubble radius computed for FC-72 was 0.05-0.1um. 


5. Effects of dissolved gases in the fluid is seen at low content values as an increase 
in run to run variability. 


6. The effect of subcooling on pool boiling incipience is apparently small. 
7, A wide variation in incipience superheat values for nominally identical runs were 
observed. 
You, Simon, Bar-Cohen and Tong [Ref. 6] conducted an experimental investigation 
of boiling in R-113. The experiment was conducted using 0.13mm diameter Chromel 
wires and 0.5I1mm diameter platinum thin film heater. The following was found. 


1. A Scanning Electron Microscope (SEM) analysis was conducted on the surface of 
the thin film heater, in which the chromel wire surface was found to contain a 
higher density of nucleation sites compared to the thin film platinum heater. 


2. The platinum heater surface was found to have a higher incipience superheat than 
the chromel wires , 73° C versus 38° C maximum. This tends to correlate well with 
item 1. 


3. Again, a wide variation in incipience superheat values for nominally identical runs 

Were experience. 

Marto and Lepere [Ref. 7] conducted experiments with surface enhancements for 
improving pool boiling performance. Marto and Lepere using FC-72 and R-113 as 
working fluids compared the performance of three commercially available boiling en- 
hancement surfaces to plain copper tubing. In each case the enhanced surface resulted 
in a lower wall surface temperature compared to the plain copper tube. However, the 
required incipience superheat for nucleation remained unaltered. 

Bergles and Kim [Ref. 2} analyzed methods to reduce the temperature overshoot of 


immersion cooling. Their method used an active bubble generator in the fluid located 


below the simulated micro chip. They used the bubble generator to activate dormant 
and extinct nucleation sites. They obtained a reduction of incipience superheat required 
from 22-33° C to 8-12° C in a working fluid of R-113. 

Maddox and Mudawar [Ref. 8] studied the effect of subcooling and surface aug- 
mentation on values of Critical Heat Flux (CHF). They found that the CHF of FC-72 
can achieve a level greater then 360 Watts/cm? with combined use of subcooling and 
surface augmentation. 

Danielson, Tousignat and Bar-Cohen [Ref. 9] conducted an analysis of a large 
number of commercially available prefluorinated inert fluids. The method of investi- 
gation used a 0.25mm diameter platinum wire heating element, with the working fluid 
at Saturation conditions. Important points reported are: 

1. The critical heat flux of FC-72 was reported as approximately 20.3 Watts/cm?. 


2. The nucleate boiling characteristic of highly wetting fluids can be predicted using 
correlations developed for water and other fluids, but are highly dependent of the 
heating surface condition. 


3. Single phase convection data are predicted accurately by a correlation developed 
by Kuehn and Goldstein.[Ref. 10] 


4. The superheat excursion at incipience boiling always appears, but the magnitude 
varies widely. 


5. The apparent radii of nucleation bubbles were calculated to be 0.1 to 0.2um. 


Kuehn and Goldstein [Ref. 10] has obtained a correlation of natural convective heat 
transfer for infinite cylinders. This correlation was used in this thesis data analysis due 
to the low rayleigh numbers obtained in the experimental runs. 

From the 3M Product Manual [Ref. 3] it can be found that the highly wetting 
dielectric fluids have a high solubility for gases. FC-72 can contain up to 48% air by 
volume compared to 1.9% air for water. This is the necessity of thoughly degassing the 


fluid prior to the experimental runs. 


D. OBJECTIVES 
The objective of this work was to investigate the effects of the boiling wake from a 
heated wire upon the heat transfer rate from wires placed within the wake. This work 


was carried out using the dielectric fluid FC-72. 


II. EXPERIMENTAL APPARATUS 


A. DESCRIPTION OF COMPONENTS 

The experimental apparatus used in this study is a redesign of an experimental ap- 
paratus used in a previous thesis. Emphasize on this design was to minimize systematic 
errors, Minimize uncertainty, and to correct deficiencies to ensure an overall increase in 
the data’s creditability over the previous apparatus, [Ref. 11]. 

In Figure 2 an overall view of the experimental setup is shown, identifying all major 
components that were utilized. Figure 3 is a photograph of the actual experimental 
chamber. 

The apparatus chamber consists of a box constructed of 0.75 inch thick 
polycarbonate sheets with a 0.25 inch thick aluminum coverplate fitted with a rubber 
O-Ring seal. Inside the box are two bulk heaters to raise the fluid to saturation tem- 
perature for degassing purposes, and to maintain the bulk fluid temperature near con- 
stant. The aluminum cover plate is mounted with three thermo-electric coolers, and 
serves as a condenser, where by the vapor condenses and returns to the chamber by 
gravity. An insert board is held in the chamber on which are mounted the test wires, 
tension springs, thermocouples, and voltage pickoffs. 

1. Test Chamber 

The interior dimensions of the chamber are 6.0 inches wide, by 2.375 deep, by 
6.25 inches high , as shown in the isometric view of Figure 8. The entire chamber was 
assembled using IPN #16 Acrylic bonding adhesive and a bead of Di-Metheyl Ethelene 
was used to seal all seams along the interior surface to ensure liquid tightness. 

To permit access for instrumentation and power leads, and as a means to ensure 
pressure equalization, a 0.75 inch diameter and 0.25 inch diameter hole were placed into 
the rear face of the chamber at the lower right and upper right corner respectively, as 
shown in Figure 4. Mounted securely around these access holes are lengths of tygon 
tubing of approximately 18 inches in length, with the opened ends suspended well above 
the chamber. 

On the interior side walls, slots are milled to permit the placement and removal 
of an instrumented wire board on which the experiments are conducted. Figure 6 is a 


photograph of the insert board, while a schematic drawing is shown in Figure 7. 


2. Aluminum Cover Plate 

A 0.25 inch thick by 7.00 inch long by 3.25 inch wide aluminum plate was uti- 
lized for use as a cover plate for the chamber, and as a condenser surface for the FC-72 
vapor. Mounted on the interior face are two thermocouples to monitor the condensor 
surface temperature, as shown in Figure 9. On the exterior face of the plate, three 
thermo electric coolers were installed for heat removal. To ensure a proper seal between 
the aluminum plate and the box structure , a shallow O-Ring groove was machined into 
the interior face and an O-Ring of 0.0624 inch diameter was installed. To adequately 
secure the cover plate to the chamber box, eight stainless steel screws of 1/16 inch di- 
ameter are used to compress the O-Ring against its mating surface. 

As mentioned previously, three solid state thermoelectric cooling devices from 
the MELCOR corp were employed to cool the cover plate condensor. They were nom- 
inally operated about the 0.4 amps and 2.0 volts operating point, from [Ref 12] their 
Operation curves and technical data is available. 

3. Heaters 

As a means of degassing the FC-72 liquid and for bulk temperature control, two 
strip heaters of 125mm x 11 mm were employed. These strip heaters were secured to the 
bottom interior surface of the chamber with epoxy, and were orientated lengthwise in the 
chamber such that one heater strip will be in front of and one heater strip behind the 
insert board. To raise the bulk fluid temperature to saturation temperature and for 
degassing purposed, both heaters are on line and operated at 1.8 amps and 20 volts. 
Once degassing is complete, the forward strip heater is electrically isolated and the rear 
heater is run at 0.9 amps and 20 volts. This will maintain the forward portion of ex- 
perimental chamber at near saturation temperature with minimum stratification of fluid, 
and maintain liquid in a degassed state since the fluid is still boiling from the heater 
strip’s surface. 

4. Thermocouples 

A total of Six copper-constantan type thermocouples with a wire diameter of 
0.05 inches (0.127 mm) were utilized for temperature measurements. Two thermocouples 
were bonded to the underside of the aluminum condensor plate utilizing OMEGA Bond 
type 101 thermally conductive epoxy, four thermocouples were installed on the insert 
board, and one was used for measuring the ambient air temperature. 

5. Insert Board 
An insert board was designed to mount into two 0.250 inch wide grooves that 


were milled into the chamber side walls. Onto this insert board are mounted four 


thermocouples to measure bulk fluid temperature, five platinum wires of nominally 0.002 
inch (0.05mm) diameter by 4.0 inch long for conducting the actual boiling experiments. 
Five springs are used to maintain the platinum wire tension so no sagging develops due 
to thermal expansion of the wires. In Appendix A the platinum wire calibrations and 


specifications are shown in detail. 


B. INSTRUMENTATION 
In the design and assembly of this experimental apparatus, five channels were wired 
and instrumentation provided to measure voltage drops at select points. Each of these 
channels had a 2.0 ohm precision resistor wired in series with the platinum wire, where 
the voltage drop across the resistor was used to obtain the current flow through the 
platinum wire. The voltage drop measured across the platinum wire is used to obtain 
the wire resistance, power dissipated and the heat flux of the platinum wire. The 
platinum wire temperature is obtained from its resistance vs temperature calibration 
curve obtained prior to experimental runs. Instrumentation was additionally used to 
obtain the bulk fluid temperature, which will be required for calculating fluid properties. 
1. Platinum Wire 
Each of the five platinum wires were run in series with a 2 ohm (+ 1%) precision 
resistor of 25 watt power rating. The following table gives the platinum wire and resistor 


combinations which were used. 


Table 1. MEASURED RESISTANCE VALUES FOR 2 OHM RESISTORS 


Wire Chan- | Measured Resistance, 20 Sample 20 Sample Std. Dev. 
nel average 


1 19980 
[a 200s ors 











A 2.0 ohm precision resistor was chosen over the previous design’s 0.1 ohms 
precision resistor in order to reduce the calculated currents uncertainty, and there by 


reduce the calculated platinum wire temperature uncertainty, [Ref. LI: p. 25]. 


2. Power Supplies 


The following power supplies were utilized and their primary function listed. 


. Hewlett Packard 6214C 0-10volts/0-1 amps: Used to power the thermal electric 


coolers. 


. Hewlett Packard 6214C 1-10volts / 0-1 amps: Used to power the constant heat flux 


platinum wire. 


. WP711 0-40 volts / 0-1 amps : Used to power the varying heat flux platinum wire. 


Hewlett Packard 6286A 0-20volts / 0-10 amps : Used to power Heater Strips for 
degassing and bulk temperature maintenance. 


Hewlett Packard 6289A 0-40 volts / 0-1.5 amps : Used to power an auxiliary illu- 
mination system or an augmenting cooling fan for the thermal electric coolers. 
3. Acquisition Unit 


To obtain the necessary measurements, a Hewlett Packard 3852A Data Acqui- 


sition System was employed, utilizing a Hewlett Packard 300 series desktop computer. 


The HP3852A utilized the following utility modules inserted into it’s chassis. 


Ir 
2. 


Hewlett Packard Model 44701A Digital Integrating Voltmeter 


Hewlett Packard Model 44705A 20 channel relay multiplexer for measuring re- 
quired voltage drops. 


. Hewlett Packard Model 44713A 24 channel relay multiplexer with electronic cold 


junction for direct conversion of thermal couple emf to temperature. 


The HP300 desktop computer was used to control the HP3852A Data acquisition sys- 


tem. Additionally, through specifically written programs, the computer can accomplish 


the required data reductions, (See Appendix D). 
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Figure 3. Experimental Chamber in Detail 
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Figure 4. 


Side View of Experimental Apparatus 
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Figure 5. Rear View of Experimental Chamber 
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Insert Board 


Figure 6. 
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Figure 7. schematic Drawing of Insert Board 
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Figure 9. Condensor Surface Details 
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Il. EXPERIMENTAL PROCEDURE 


A. PREPARATIONS 
Careful preparations were required to ensure that the data obtained was of high 
quality. The most notable requirement being the necessity of thoughly degassing the 
FC-72 liquid, [Ref, 3: p. 3]. 
1. Normal Warm-Up Procedure 
In carrying out the experimental runs the following procedures were followed. 


1. Power-up the HP3852A Data Acquisition Unit to start the warm-up period of one 
hour minimum. [Ref. 13: p. 586] 


2. Ensure the liquid level is above Wire 5, but below the condensate drip tray. 
3. Set power supply of thermoelectric coolers to 1.3 amps and 2.0 volts. 


4. Ensure the experimental chamber’s insulation jacket is in place and set power sup- 
ply to strip heaters to approximately 1.8 amps and 20 volts. 


5. Verify the 0.25 inch diameter vent line is unobstructed. 


6. Power on the HP300 series desktop computer, and load data acquisition program 
as contained in Appendix D. 


7. Once fluid has reached saturation and is boiling vigourisly, hold for a minimum of 
one hour, typically two hours was achieved, to permit adequate degassing of fluid. 
In actual practice, the time required to raise the bulk temperature of the fluid to satu- 
ration temperature of 56 °C was nominally two hours. 
2. Post Warm-up Procedure 
Once the required degassing has been accomplished, power is temporarily se- 
cured to the strip heaters. Once power is secured the forward strip heater is electrically 
isolated. The rear strip heater is powered at 0.9 amps and 20 volts, and the thermo 
electric coolers are secured. 
The power setting of the strip heaters permits vigorous boiling from it’s surface. 
This serves to keep the fluid degassed and the bulk fluid temperature near saturation and 
near isothermal. By keeping the fluid level below the insert board’s upper edge, the fluid 
forward of the insert board remains undisturbed by the boiling action occurring in the 


rear chamber. 
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DATA ACQUISITION PROCEDURE 


Upon completion of the previous preparatory procedures, the experimental runs are 


ready to begin. The procedure followed for an experimental run is as follows: 


l. 


Plug power supply into power receptacles for wire 1. Power on and ensure current 
limiter is set to max value of 1.0 amps. 


. Plug HP6214C power supply into power receptacles of desired constant heat flux 


wire. 


3. Set power supply of step 2 to proper output to give desired heat flux from the wire. 


tf 


c. 


. Start running the Data Acquisition Program, which is programmed to measure the 


following data from the desired wires at two minute intervals. 
e Precision Resistor Voltage drop, up to five channels. 

e Platinum Wire Voltage drop, up to five channels. 

e Four bulk fluid temperature measurements 

e Two condenser surface temperature measurements. 

e One ambient air temperature measurement. 


After each data interval, increase the voltage to the varying heat flux wire by 0.2 
Volt increments. This was continued until the on line power supply was at 2.3 volts 
of outputs, the increments were then reduced to 0.1 volts until nucleation occurred. 
The increments were then returned to 0.2 volts until a heat flux of approximately 
200,000 watts/m? were achieved. At any higher heat flux the wire will go into the 
film boiling regime. 


. Reduce the power to the varying heat flux wire at 0.3 Volts increment to bring the 


wire back to a no power condition. 


Repeat steps 3 thru 6 for all other wire and heat flux configurations. 


DATA REDUCTION 


Data that was obtained from the actual experimental apparatus were two voltage 


measurements per active platinum wire channel, and six temperature measurements. 


The process used for data reduction is a combination of computer programs and 


graphical outputs for numerical presentation and visual presentations respectively. 


l. 


Z 


1. Numerical Data 


The numerical data values required to be processed are: 
Fluid Bulk Temperature 


(7, eee eae 
Tue = ag - - (-G) (3.1) 


and where, T, thru T, are the thermocouple measurements of the bulk fluid. 


Current flow through the Platinum Wire, where 
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V 
I “<5. (Amps) (3.2) 


wire — R 





22 
3. Platinum Wire's electrical resistance value, where 
V prwire 
Re:wire = —___ (Ohms) (3.3) 
wire 


4. Platinum Wire surface temperature, where 
T sur = ER pr wire aE (Age) (3.4) 
where, 
T,,,= the actual wires surface temperature 


€ is the slope of the calibration curve 
T, 1s the calibration curves y-axis intercept 


5. Fluid Film Temperature, where 


Toe + Upeik 
T tm an i u (3.5) 
6. Dissipated Heat Flux of platinum wire, where 
Lite x V prwire Watts 
= a el ee ae (3.6) 


Alea gyrr m 


where, Areag,,,= 7 X length,,,. X Dur 


The above numerical data were all derived from the raw data, and previously 
determined calibration curves. The following table provides physical dimensions and 


temperature calibration coefficients of the platinum wires. 


Table 2. PLATINUM WIRE PHYSICAL PARAMETERS FOR EQUATION 
(3.4) AND EQUATION (3.5) 


[Wie] | + enh —«~Sdiamerer 
oe | co | ce | mm | om 
[7302s 9922 —< os 
[4 3.asa_ | asses. =| ios 









To see in detail the methodology used to obtain the numerical values presented 


in the preceding table, proceed to Appendix A ‘Platinum Wire Calibration’. 
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2. Fluid Property Data 
The Fluid Properties to be determined are as follows: 


. Thermal Conductivity of fluid, as obtained from Figure 5 of the 3M Corporation 
Flourinert Product Manual [Ref. 3: p. 15] have been determined to be, 


0.6033 — 0.00115 x T : 
flim , Watts 
ia | 10 eel mOnL Fe 


. Thermal Expansion Coefficient, from table 4B in the Product Manual [Ref. 3: p. 10] 
is given as 
0.00261 
= 1.740 — 0.00251 x Thin (<E = ) oY 
. Kinematic Viscosity of the Fluid , as obtained from Figure (3) in the Product Man- 
ual [Ref 3: p.13] and determining an exponential curve fit yields 


K m? 


v = 1.203952E —8 x exp(1,058.4109 T+ 2731sc hs? 


(3.9) 


. Specific Heat of the Fluid for all the Flourinert Fluids are obtained from Figure (4) 
of the Product Manual [Ref 3: p. 14] 1s, 


= (0.24111 + 3.70374E — 4 X Tyan) x 4186 a (3.10) 


. Liquid Density of FC-72, from expression in table 4B and constants from table 4C 
in the Product Manual [Ref 3: p 10] 1s given by 





k 
p = (1.740 — 0.00261 Ty) x 1000 —=) (3.11) 
m 
. Thermal Diffusivity of the liquid is given by, 
maak 
o= px, (3.12) 
. Prandlt Number of the Fluid is given by 
Pr=—— (3.13) 
. Ralyeigh Number of Fluid based on wire diameter is given by, 
x Bx (T,,-— T x D? 
Rae g p ( sur ol) (3.14) 


Vx a 


. Correlation for the Nusselt Number was obtained by using the correlation of Kuehn 
and Goldstein [Ref 10: p. 1128] given by 


ja 
2 


eee 
[(0.518 x Ray x [1 + (S24 P95 + 0.1 x Ray?) Fy"? 


Nup = 
Inf] + 
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10. Theoretical Heat Flux as obtained from using the N, value obtained from equation 
(3.15). 


k x Nup A ee eats 
gq = RNY * Tou = Tout) ¢ Watts 


- (3.16) 


m 

With all the required numerical data and fluid properties obtained, it is now 
necessary to present the numerical information in a graphical representation to see the 
relevancy of the results obtained. This will become apparent with the presentation of 
results in Chapter IV ‘Results and Discussion’. 
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IV. RESULTS AND DISCUSSION 


The data obtained can be grouped for discussion as follows: 
1. Individually powered wires 


Saturated Conditions 
Sub-Cooled Conditions 


2. Bubble Pumping Effect Runs 


Saturated Conditions 
Sub-Cooled Conditions 


A. INDIVIDUALLY POWERED WIRES 

Before analyzing the pumping effect of the bubbles on the upper wires, the individ- 
ual wire boiling experiments will be investigated. All of the individual powered wire runs 
shown here were made at or near saturated fluid conditions with two exceptions. The 


following table will provide a listing of single wire runs. 


Table 3. INDIVIDUAL WIRE RUNS 


ae 
3 2 















In all cases conditions were under atmospheric pressure with 0.05mm diameter platinum 
wires. 
1. Natural Convection 
Figure 10 is a typical boiling curve for FC-72 with the Kuehn and Goldstein 
natural convection region shown along with the experimental points. The nucleate 
boiling zone starts at boiling incipience, which is evident from the sudden drop of wall 


temperature, a characteristic of a highly wetting fluid. Figure 11’is a compilation of wire 
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] data, at or near saturation, obtained over a period of time. As can be seen, with the 
exception of the boiling incipience point, the data is virtually super imposable over each 
other. This reproducability of data over time indicates no significant aging has occured 
on wire I, and more importantly no change in calibration. 

In fact in the natural convection region Figure 10 thru Figure 16 are all essen- 
tially superimposable over each other. Figure 17 is an overlay of all the previous fig- 
ures. Therefore, it appears that whether the fluid is Saturated or Sub-Cooled, the fluid 
bulk temperature appears to have little effect on the natural convection heat transfer 
rate. 

There was generally very good agreement between the natural convection data 
and the Kuehn and Goldstein correlation [Ref 10] (equation 3.15). The only deviation 
existing between the figures being the total temperature excursion to incipience, which 
ranged from 27° C to 38° C. 

2. Nucleate Boiling 

As previously stated, the nucleate boiling zone starts at boiling incipience. Us- 
ing wire | as a standard and Figure 10as base line boiling characteristic data a compar- 
ison to other wires at saturation condition, and a comparison to subcooling effects on 
wire’s 1 boiling can be made. As is evident from Figure 17, the boiling characteristic 
of all the wires at saturation condition tend to overlap each other. The differences aris- 
ing tending to be caused by individual wire placement in the fluid. However, the effect 
of subcooling does result in a decrease in the boiling curve overshoot at point of 
incipience. This decrease in the boiling curve overshoot results in a net decrease in the 
superheat required to initiate nucleate boiling. 

3. Departure from Nucleate Boiling 

The Critical Heat Flux (CHF) for FC-72 at saturated condition was calculated 
by Danielson et all [Ref 9] to be approximately 203,000 Watts/m?. Experimental value 
for CHF ranged from 198,000 Watts/m? to 220,000 Watts/m’, a spread of -2.5% to 
8.4%. A CHF for a subcooled condition of T,,,= 38° C was found experimentally to 
be approximately 400,000 watts/m?. A CHF for T,,,= 23 ° C was not obtained during 


the experimental runs. 
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HEAT FLUX (WATTS/M2) 


Figure 10. 


WIRE 1 BOILING CURVE 


NEAR SATURATION T,,, = 55.5° Celisus 


LEGEND 
0 = HEAT FLUX1 INCREASE 
o = HEAT FLUX1 DECREASE 
GOLDSTEIN CORRELATION 


20 Zs 30 
TSURF — TSAT (C) 


Wirel boiling curve T,,, = 55.5°C 
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WIRE 1 BOILING CURVES 


ALL BOILING CURVES AT SATURATION OR NEAR SATURATION 


100000 


HEAT FLUX (WATTS/M2) 


LEGEND 
a = WIRE 1 DATE 14 AUG 
o = WIRE 1 DATE 20 AUG 
a = WIRE 1 DATE 19 AUG 
GOLDSTEIN CORRELATION 
1000 
25 30 35 40 
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TSURF — TSAT (C) 


Figure 11. Wire 1 Boiling Curves over time 
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WIRE 1 BOILING CURVE 
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Figure 12. Wire 1 boiling curve T,,, = 38°C 
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Figure 13. 


WIRE 1 BOIUNG CURVE 


DUBGOOLED 7. 7— 22.5 Celsus 


LEGEND 
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GOLDSTEIN CORRELATION 
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Wire | Boiling Curve T,,,, = 23°C 
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WIRE 2 BOILING CURVE 


SATURATED Dux = 56.0° Celisus 


HEAT FLUX (WATTS/M2) 


LEGEND 
O = HEAT FLUX2 INCREASE 
o = HEAT FLUX2 DECREASE 
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Figure 14. Wire 2 boiling curve T,,, = 56°C 
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Figure 15. 


WIRE 3 BOILING CURVE 


SATURATED T,,,, = 56.0° Celisus 
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WIRE 4 BOILING CURVE 


SATURATED T,,,, = 56.0° Celisus 


HEAT FLUX (WATTS/M2) 
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Figure 16. Wire 4 boiling curve T,,, = 56°C 
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Figure 17. 
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B. BUBBLE PUMPING EFFECT ON HEAT TRANSFER RATE 

In the upcoming discussion a brief explanation of the figures are required. In each 
case for these experiments, heater wire number | was used to generate the bubbles for 
examining the pumping effects upon the upper wires. The figure insert is used to show 
which wires are being utilized. Only select data points were plotted for the passive 
(constant heat flux wire) wire, and those points correspond in time to the identically 
numbered data point on wire 1’s boiling curve. The data points were chosen since they 


define regions of wire 1’s boiling curve as shown in the following table. 


Table 4. DATA POINT DESCRIPTION 
| Natural Convection Start 
po 8 Partin! Boiling om wire 
po Full nucleate Boiling 
pS Max Nucleate Boiling 
pS End often 









1. Passive Wire at Low Heat Flux 
In these experimental runs the heat flux on the passive wire was set to approx- 
imately 8500 watts/m?. At this low heat flux it was expected that little to no active 
nucleation will occur on the passive wire, even after wire 1 is in nucleate boiling. 
Therefore, all the enhancement of the upper wire will be due to an addition of a force 
convection term to the heat transfer rate. 
In comparing Figure 18 to Figure 19the following can be inferred: 


1. Total reduction of wall temperature on the passive wires are comparable. Wire 2 
having a 3.7°C and Wire 4 having a 3.0° C drop in wall temperature. 


2. A definite effective distance exist over which wire’s 1 influence 1s imparted upon the 
upper wires, for example: 


Wire 2 and wire 4 showed a very small change in wall temperature at point 2 (point 
of incipience). 

After incipient boiling has occurred on wire 1, wire 2 showed a greater drop in wall 
surface temperature then wire 4. Approx. 2.5° C drop for wire 2 and 0.1° C drop 
for wire 4, 


by. 


3. At point 4, the point of maximum temperature in nucleate boiling, both wire 2 and 
wire 4 showed a comparable surface temperature, approximately 2° C above satu- 
ration. 

As can be seen from the figures, some nucleation sites have been activated on wire 2 and 
wire 4, with point 5 showing a lower wall temperature then point 1. 
2. Passive Wire at Medium Heat Flux 

In these experimental runs the heat flux on the passive wire was set to approx- 
imately 20,500 watts/m?. At this heat flux it was expected that active nucleation will 
occur on the passive wire after wire | is nucleate in boiling. 

In comparing Figure 20 to Figure 21 
the following can be inferred: 


1. The total reduction of wall temperature on the passive wire , points 1 to points 4, 
are comparable at 7.1° C for wire 2, and 6.3° C for wire 4. 


2. The spatial effect is evident on the influence of wire 1 on the upper wires, for ex- 
ample: 


Between point 1 and point 2, wire 2 shows a greater reduction in wall surface tem- 
perature than is evident in wire 4. 


At point 3, wire 2 shows a much greater drop in surface temperature than wire 4: 
4° C compared to 2° C. The difference being the amount of 


3. At point 4, maximum nucleate boiling from wire 1, both wire 2 and wire 4 show 
comparable values for surface temperature, approximately 
les =... 
As can be seen from the figures, wire 1 has activated previously inactive nucleation sites 
on wire 2 and wire 4. This is evident from point 5 being at a lower surface temperature 
than point 3. 
3. Passive Wire at High Heat Flux 
In these experimental runs the heat flux on the passive wire was set to 40,000 
Watts/m? on wire2 and 35,000 Watts/m? on Wire 4. At this heat flux it was expected 
that a large number of inactive nucleation sites will be activated after wire | is nucleate 
in boiling. 
In comparing Figure 22 and Figure 23 it is evident that the information ob- 
tainable from the figures is comparable to the information for the medium heat flux case. 
However, the wall temperature reductions on the passive wires are much greater as 


compared to the Medium Heat Flux case, in particular the following: 


1. In comparing the surface temperature reduction from point 1 to point 2, wire 2 has 
dropped 2.5° C while wire 4 had almost a zero drop of temperature. This appears 


a, 


to be due to the relative locations of the passive wires with respect to wire 1, and 
is especially apparent at the higher heat fluxes. 


2. Many inactive nucleation sites have been activated by wire 1’s nucleation bubbles. 
As can be seen by comparing point | to point 5(6), a difference on Wire 2 of 13° 
C and on Wire 4 of 7° C is apparent. 
4. Wire 1 Passive with Wire 2 varying 
In these experimental runs, wire | was set to 150,000 Watts/m? so it was boiling 
vigorously. wire 2 was actively varied over a large heat flux range as can be seen from 
Figure 24. On Figure 24 a previously obtained boiling curve of Wire 2 was overlayed 
onto it. Upon comparing the two curves the most striking feature is the disappearance 
of the Hysterises Loop as is evident in the overlayed boiling curve data. In this ar- 
rangement, the boiling curve overshoot has been removed as a design problem. bubbles 
from Wire 1 has removed a major design consideration. 
5. Subcooled 
In these experiments, the bulk fluid was degassed and then subcooled to a tem- 
perature of 38° C. In these experimental runs, wire 2 was set at 35,000 watts/m? and 
53,000 watts/m?, as shown in Figure 25 and Figure 26 respectively. In subcooling, no 
significant effects were notable beyond those seen in the saturated runs. The only effect 
that was apparent was the ability to run wire 2 at a greater heat flux with out incipent 


boiling occuring on the wire. 
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WIRE 2 AT 8,500 WATT/M2 


NEA oS a URATED [7 — 55,5 Celsius 


HEAT FLUX (WATTS/M2) 


LEGEND 
“10 = HEAT FLUX1 INCREASE 


o = HEAT FLUX2 CONSTANT 
6 = HEAT FLUX1 DECREASE 
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Figure 18. Wire2 at 8500 W/ m? T,,, = 55.5°C 
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WIRE 4 AT 8,500 WATT/M2 


NEAR SATURATION T,,y, = 55.6° Celsisus 


HEAT FLUX (WATTS/M2) 
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Figure 19. Wired at 8500 W/ m T,,,,, = 55.6°C 
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HEAT FLUX (WATTS/M2) 


Figure 20. 





WIRE 2 AT 20,500 WATT/M2 


NEAR SATURATION T,,,, = 55.6°Celisus 
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Wire 2 at 20,500 W/ m? T,,, = 55.9degreeC 
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WIRE 4 AT 20,500 WATT/M2 


NEAR SATURATION T,.,, = 55.9° Celisus 
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Figure 21. Wire 4 at 20,500 W/ m? T,,, = 55.6°C 
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WIRE 2 AT 40,000 WATT/M2 


NEAR SATURATION T,,, = 54.6° Celisus 


HEAT FLUX (WATTS/M2) 


LEGEND 
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Oo = HEAT FLUX2 CONSTANT 
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Figure 22. Wire 2 at 40,000 W/ m: Tye = 94.6°C 
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WIRE 4 AT 35,000 WATT/M2 


NEAR SATURATION T,., = 55.0°Celisus 
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Figure 23. Wire 4 at 35,000 W/ m? T,,., = 55°C 
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WIRE 1@ 180,000 WATTS/M2 


SATURATED 7,,,=56°C 
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Figure 24. 
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Wire 1 passive, Wire 2 active T,,, = 56°C 


41 





35 


40 


HEAT FLUX (WATTS/M2) 


Figure 25. 


WIRE 2 AT 35,000 WATT/M2 


SUBCOOLED 7,.,, = 38.0° Celsisus 
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Wire 2 at 35,000 W/ m T,,,, = 38°C 
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HEAT FLUX (WATTS/M2) 


Figure 26. 


WIRE 2 AT 53,000 WATT/M2 


SUBCOOLED T,,, = 38° Celisus 
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Wire 2 at 53,000 W/ m? T,,,, = 38°C 


43 





33 


40 


V. CONCLUSIONS 


The study conducted in this thesis was to investigate the effects of the wake from a 


bubble generating wire on wires located downstream within the wake. Experimental 


runs were conducted at saturated condition and subcooled conditions in which the fol- 


lowing were found: 


A. SINGLE WIRE BOILING 


l. 


The effects of subcooling on the natural convective heat transfer from the wire was 
negligible. 


The effects of subcooling tended to increase the heat flux required for incipient 
boiling to initiate. 

With subcooling, boiling hystersis is still observed, but the hystersis magnitude de- 
creases as subcooling increases. 


Subcooling results in a reduction in wall temperature after the start of nucleate 
boiling as compared to saturated cases. 


B. PASSIVE WIRE ANALYSIS 


Heat Transfer enhancement effects from Wire 1 upon the downstream wires has 


many factors within it. The following factors have been identified: 


l. 


Ih. 


An enhancement to convection resulting from buoyancy induced fluid flow off of 
wire 1 was seen to exist. This effect is very spatially dependent with the effects 
most pronounced on wire 2 and negligible on wire 4. This corresponds to point 1 
and point 2 on the figures. 


The sudden drop of wall temperature between point 2 and point 3/4 was seen to 
be the results of the following: 


Nucleation bubbles from wire | activating the inactive nucleation sites on the upper 
wires resulting in boiling heat transfer off of the upper wires. 


An additional force convection enhancement due to the nucleation bubbles from 
wire | breaking up the thermal boundary layer about the upper wires is believed to 
exist. 


The drop in temperature on the upper wires corresponding with points 3/4 to point 
4/5 appears to be due purely to an increase in force convection due to the high 
density of bubbles off of wire 1. 


Additional conclusions reached are. 


A constant source of bubble generation from below is an effective way to eliminate 
the Boiling Curve Overshoot and Hysteresis loop associated with dielectric fluids. 


2. At low heat fluxes, proximity to Wire | is important for heat transfer enhancement 
to occur. This dependency reduces as heat flux increase on the passive wire. 
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APPENDIX A. CALIBRATION 


A. PLATINUM WIRE 


1. Surface Micro-Structure 


Prior to the experiments, a length of 0.05mm platinum wire was analyzed using 


a Scanning Electron Microscope (SEM). The length of wire analyzed was not used in 
any experiments, but it did come from the same spool of wire as the experimental lengths 


of wire, therefore, it should be a representative sample. In analyzing this sample the 


SEM brought out several important features: 


I. 


From Figure 27 a wire diameter of 0.052mm can be determined. Also, the surface 
is relatively unmarked at this magnification of 1490x; little scoring on the wire and 
a low density of surface pits are evident. 


In Figure 28, an increase in magnification to 4740x within the field of Figure 27, 
a more detailed view of the larger surface defects are evident. Maximum defect size 
is approx 2um in length and 0.5um in width. A surface defect of that size would 
be too large to be a possible nucleation site , according to the calculations done in 
Ref. 4 of 0. Lym. 


In Figure 29, a further increase in magnification to 13,200x within the field of 
Figure 27, pits of approximately 0.8um max down to 0.3um are seen. It will be 
these surface defects that boiling will initiate from. The density of these pits along 
the length of the wire analyzed was not great, this implies that the incipient super- 
heat will be large. 





Figure 27. Surface micro geometry: 
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Figure 28. Surface micro geometry: Sample Platinum Wire (4740x) 





Figure 29. Surface micro geometry: Sample Platinum Wire (13,200x) 
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2. Desirability for Using Platinum 
Platinum was chosen for several reasons: 
1. It has a highly linear response of Temperature vs Resistance. 


2. It has a large AR for a given AT. 


In each aspect, platinum performed very well. As can be seen in the calibration curves 
of Figure 31 thru Figure 40 response is linear, with a 0.1Q change in resistance for ev- 


ery 5° C change in temperature. 


B. BACKGROUND 
An important step in minimizing experimental data uncertainty is an accurate cali- 
bration of the platinum wires used. In obtaining the calibration curves for the wires, 


three calibration runs were performed. 


1. An initial calibration in a heated calibration bath with a platinum resistance 
thermometer. 


2. An in place calibration of wires in the test chamber, with the four installed 
thermocouples used for temperature measurements, was conducted to ensure the 
calibration curves from the calibration bath were not altered from handling the wire 
insert board. 


3. A final in place calibration of wires in the test chamber, 
with the four installed thermocouples used for temperature measurements, was 
conducted to ensure the calibration curves used for data acquisition were not al- 
tered due to any reason. 


The second calibration run was required to ensure no distortion occurred to the wires 
when the insert board, on which the wires are mounted, was removed from the calibrated 


bath and installed into the experimental test chamber. 


C. INITIAL CALIBRATION 


The initial calibration was conducted using the following equipment. 


1. Rosemount Engineering Co. Model 913A calibration bath with Etheylene Glycol 
as the working fluid. 


2. Rosemount Engineering Co. Model 923B power supply, controlling calibration 
bath temperature. 


3. Rosemount Engineering Co. Model 920A commutating bridge, used to measure an 
resistance to 0,00]Q accuracy. 


4. Rosemount Engineering Co. Model 162C S/N 985 Platinum Temperature Probe. 
Used to measure calibration bath temperature through the commutating bridge. 


5. Hewlett Packard Data Acquisition System. Used to measure platinum wire resist- 
ance using a 4-wire ohm measurement technique. 
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1. Calibration Heating Bath 
The present investigation was conducted using an externally heated calibration 
bath with etheylene glycol as the working fluid. This fluid was circulated by a centrifugal 
pump to ensure a uniform bulk temperature throughout the bath. Item (1) and item (2) 
are utilized . 
2. Temperature Measurements 
Temperature measurements for calibration were conducted using a high sensi- 
tivity commutating bridge to measure the resistance of a platinum resistance temper- 
ature probe. This platinum temperature probe had been previously calibrated to a 
standard that is traceable to the National Bureau of Standards. 
3. Resistance Measurements 
In calibrating the platinum wire the resistance of the wire was measured with the 
same Data Acquisition Unit Components as would be used in the experimental runs. 
This was to ensure the uncertainty of measured resistance would be the same as between 
the calibration run and experimental runs. 
The Data Acquisition Unit Component used are: 
HP3852A Control Unit 
HP44701A Integrating Voltmeter card 
HP44705A 20 channel Multiplexer card 
HP300 series desk top computer 


ee ae oe 


HP44713A 24 channel multiplexer with electronic cold junction for thermo couple 
measurements. 


4. Calibration Procedure 
By using the previously described components for measuring temperature and 


resistance, the following procedure was used to obtain the system’s calibration data: 
I. Set calibration bath controller for 35° C 
2. Allow bath temperature to stabilize. 


3. Carefully place Insert Board into bath liquid. EXTREME care must be taken to 
ensure that none of the wires are deformed or destroyed. 


4. Allow insert board to come to thermal equilibrium with the bath, 
this takes approximately five minutes. 


5. Take temperature reading with platinum wire temperature probe and commutator 
bridge, and resistance measurement of wires with data acquisition unit. 


6. Remove insert board, and set bath temperature controller for a 5° C increase. 


7. Repeat step 2 thru 6 until an 90° C reading has been achieved. 
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8. Activate External Cooling system. 
9. With temperature control setting decreasing at 5° C intervals, repeat step 2 thru 6 
until a 35° C reading has been achieved. 
D. IN PLACE CALIBRATION 
To ensure that the calibration curve of platinum wire | thru 5 are reliable, an in 
place calibration run was performed on the wires. This calibration procedure was per- 
formed before any experimental runs were conducted, and after the last experimental run 
was completed to detect any changes in calibration due to wire deformation or aging. 
The equipment utilized for the inplace calibration was the same Data Acquisition 
equipment used for the resistance measurements of the Bath Calibration runs with two 
additions. 
1. HP44713A 24 Channel Multiplexer with electronic cold junction. 


2. Four type T thermocouples. 


The thermocouples are now utilized, instead of the platinum resistance temperature 
probe, for measuring fluid bulk temperature. The uncertainity associated with the 
thermocouples is 0.4°C as compared to the platinum resistance temperature probe used 
as a standard. Figure 30 shows the thermocouples uncertainity compared to the 
platinum resistance temperature probe. 
1. Inplace Calibration Procedure. 
The following procedure is used for inplace calibration of the platinum wires: 
1. Ensure Data Acquisition system is wired for 4-wire ohm measurements. 
2. Turn on Data Acquisition system and allow one hour of warm-up time. 


3. Load Calibration Acquisition Program into the HP-300 series Computer. See Ap- 
pendix D for program listing. 


4. Set power supply to bulk heaters to 20volts and 1.8 amps. 


5. Start Acquisition program. The acquisition program will measure the resistance 
of each platinum wire, and measure the temperature output of each thermocouple 
at a 10 minute interval. 


6. Allow acquisition program to run until bulk fluid has reached Saturation temper- 
ature of 56° C, 
FE. CALIBRATION RESULTS 
1. Platinum Wire | 
From Figure 31 the importance of conducting an in place calibration prior to 


any experimental runs becomes obvious. wire |’s Bath calibration data has been con- 
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siderably altered while it was being inserted into the test chamber. The initial in place 
calibration data was used for the experimental runs calibration curve. The resulting 


calibration was: 


T, = §3.902 x R, a 254.470 (eC) 


From Figure 32 it can be seen that the post calibration curve matches well with 
the data acquisition curve. 
2. Platinum Wire 2 
From Figure 33, it is apparent that the initial in place calibration curve fairs in 
well with the calibration bath’s decreasing data. The initial in place calibration data was 


used for the experimental runs calibration curve. The resulting calibration was: 


T> = 53.666 x R, — 253.009 (°C) 


From Figure 34 it can be seen that the Post run calibration matches well with 
the Data Acquisition Calibration curve. 
3. Platinum Wire 3 
From Figure 35 it is apparent that wire 3’s Bath calibration data has been al- 
tered between the calibration bath and insertion into the test chamber. The initial in 
place calibration data was used for the experimental runs calibration curve. The result- 


ing calibration was: 
T; = 53.141 x R3 — 252.729 (°C) 


From Figure 36 it can be seen that the Post run calibration matches well with 
the Data Acquisition Calibration curve. 
4. Platinum Wire 4 
From Figure 37, it is apparent that the initial in place calibration curve fairs in 
well with the calibration bath’s decreasing data. The initial in place calibration data was 


used for the experimental runs calibration curve. The resulting calibration was: 


T, = 53.484 x R, — 252.848 (°C) 


Form Figure 38 it can be seen that the post run calibration matches well with 


the data acquisition calibration curve. 


a1 


5. Platinum Wire 5 
From Figure 39 it is clear to see the shift in calibration curves between the cal- 
ibration bath data and initial inplace calibration data. The calibration data has been 
altered while it was being inserted into the test chamber obviously. The initial data 
calibration was used for the experimental run calibration curve. The resulting cali- 


bration was: 


T; = 53.085 x Rs — 260.049 (°C) 


From Figure 40 it can be further seen that a considerable error has been intro- 
duced into wire 5’s calibration, an error of 3.7° C. This places an uncertainity value onto 


wire’s 5 temperature data of at least 3.7° C, an unacceptable value. 
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Figure 31. Calibration of Platinum Wire 1: Calibration curve of wire 1 of Cali- 
bration Bath Temperature Increasing, Calibration Bath Temperature 


Decreasing, and of the Initial in place calibration data. 
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Figure 32. Calibration of Platinum Wire 1: Calibration curve of wire | of Initial 


in place calibration curve, and of Final in place calibration curve 
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Figure 33. Calibration of Platinum Wire 2: Calibration curve of wire 2 of Cali- 
bration Bath Temperature Increasing, Calibration Bath Temperature 
Decreasing, and of the [Initial in place calibration data. 
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Figure 34. Calibration of Platinum Wire 2: Calibration curve of wire 2 of Initial 


in place calibration curve, and of Final in place calibration curve 
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Figure 35. Calibration of Platinum Wire 3: Calibration curve of wire 3 of Cali- 
bration Bath Temperature Increasing, Calibration Bath Temperature 


Decreasing, and of the Initial in place calibration data. 
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Figure 36. Calibration of Platinum Wire 3: Calibration curve of wire 3 of Initial 


in place calibration curve, and of Final in place calibration curve 
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Figure 37. Calibration of Platinum Wire 4: Calibration curve of wire 4 of Cali- 
bration Bath Temperature Increasing, Calibration Bath Temperature 


Decreasing, and of the Initial in place calibration data. 
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Figure 38. Calibration of Platinum Wire 4: Calibration curve of wire 4 of Initial 


in place calibration curve, and of Final in place calibration curve 
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Figure 39, Calibration of Platinum Wire 5: Calibration curve of wire 5 of Cali- 


bration Bath Temperature Increasing, Calibration Bath Temperature 
Decreasing, and of the Initial in place calibration data. 
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Figure 40. Calibration of Platinum Wire 5: Calibration curve of wire 5 of Initial 


in place calibration curve, and of Final in place calibration curve 


F. SUMMARY 
The following points were brought out during the calibration process. 


1. The process of moving the platinum wires from the calibration bath and installing 
them into the test chamber resulted in a shift in the calibration curves. 
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2. The most probable cause of all the calibration shifts were the lead wires in the back 
of the insert board being pulled and adding an additional tension onto the platinum 
wires, thereby stretching them. The fact that the calibration curves consistently 
shifted to an increase in resistance supports this theory. 


3. The maximum difference observed between the Platinum Wire Temperature Probe 
and the Type T thermocouples were 0.4° C. Therefore, a temperature uncertainity 
for the in place calibration will be set at 0.4° C. 


4. The differences between the wire’s data acquisition calibration curve and the wire’s 
final calibration curve did not exceed the thermocouples uncertainity in temper- 
ature. 


APPENDIX B. SAMPLE CALCULATIONS 
In these sample calculations the output value of wire 1 will be used. 


A. DETERMINATION OF DERIVED VALUES 
1. Determination of the Average Bulk Temperature 
The average bulk temperature of the liquid was obtained by the arithmetic av- 


erage of the four immersed themocouple outputs. 


(ay, + T, +T; + T,) 


T pute aa 4.0 
5 (55.835 + 55.884 + 55.988 + 55.932) 
a 4.0 
= 55.9098°C 


2. Determination of Wire Current 
The current flow through the platinum wire is obtained by measuring the volt- 


age drop across a precision 2 ohm resistor. 


I = V Res 
PtWire — 
Ry 


_ (0.74605V olts) 
~ (1.99940hm) 
= 0.2381Ohms 


B. DETERMINATION OF PLATINUM WIRE RESISTANCE 
the resistance value of the platinum wires are obtained by measuring the voltage 


drop across the desired wire and dividing it by the current flow through it. 


Vp, Wire 


Rorywire a I 
PtWire 


— (1.7233Volts) 
~ (0.2831Amps) 
= 6.0874 Ohms 


1. Determination of Wire Surface Temperature 
The wire’s surface temperature are obtained from the calibration curves shown 
in APPENDIX A, wire | will be shown as an example. 
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T prwiret == 53.902 x Rp wire — 254.470 


= 53.902 x 6.0874 — 254.470 
= 73.356°C 


2. Determination of Input Power 
The input power of the wires are determined from the current flow through the 


wire times the voltage drop across the wire. 


Power, = Vppwirer X Uppyires 


= 1,.7233Volts x 0.28310Ohms 
= 0.48787 Watts 


3. Determination of Heat Flux 
The heat flux dissipated by the platinum wire is determined from the inputted 


power divided by the wire’s surface area. 


Power, 
A, 
(0.48787 Watts) 
1.5585E-5m? 
= 31,303.44Watts/m” 


Hflux, ae 





where, 


A,=2zxD,xL, 
= 7X 5E-5 m x 0.09922m 
= 1,5585E-5’m? 


C. DETERMINATION OF FLUID PROPERTIES 
1. Film Temperature 
The film temperature is the arithmetic average of bulk fluid temperature, and the 


Wire's surface temperature. 


T _ (Tyyiea! pe Wire) 
fim 2.0 


(54.0°C + 58.0°C) 


2.0 
= 56.0°C 


66 


2. Thermal Conductivity 
From Figure 5 of the 3M Corporation Flourinert Product Manual [Ref. 3: p. 


15], the thermal conductivity has been determined to be: 


_ (0.6033 — 0.00115 x Tim) 


10. 
= 0.05389 Watts/m°C 


at Trim = 56°C 
3. Liquid Density 
Using the expressions in Table 4B and the constants presented in Table 4C of 
the 3M Corporation Product Manual,[Ref. 3: p.10], the liquid density has been deter- 


mined to be: 


p, = (1.740 —0.00261 x Tym) x 1000 
= 1,593.8kg/m° 


at Tym = 56°C 
4. Kinematic Viscosity 
Using Figure 3 in the 3M Corporation Product Manual, [Ref. 3: p.13] and de- 


termining an expodintial curve fit valid from 0°C to 90°C yields: 


= oe eee 
2S Pe TS. 


A = 1.203952E — 8m’/s 
B = 1,058.4109°K 


at Trim = 56°C 
v = 3.000E — 7m?/s 


5. Specific Heat 
From Figure 4 of the 3M Corporation Flourinert Product Manual, [Ref. 3: 
p.14], the specific heat has been determined to be: 


Cp = (0.241111 + 3.70374E — 4 x Tym) x 4186 
= 1,096.15 /kg°C 


at Tam = 56°C 
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6. Thermal Expansion Coefficient 
Using the expression in Table 4B and the constants presented in Table 4C of the 
3M Corporation Product Manual, [Ref. 3: p.10], the thermal expansion coefficient has 


been determined to be: 


B= 0.00261 
(1.740 —0.00261 x Tym) 


= 0.00164 1/°C 


at Trim = 56°C 
7. Thermal Diffusivity 
Is determined to be at Ty, = 56°C 


2 Sse 
px Cp 


= 3.0848E — 8m’/s 


8. Prandlt Number 
Prandlt Number is determined to be at Ty, = 56°C 


|< 


|e 


I 
Sg 
~~] 
RQ 
Gn 
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APPENDIX C. UNCERTAINTY ANALYSIS 


A. UNCERTAINTY IN SURFACE AREA 


Wn = 9.003 mm 
Oo, = 0.8 mm 
A = nxDL 
where: 
D = Wire Diameter 
L = _ Length of platinum wire 
A = _ Surface Area of platinum wire 
OA 
apie oe 
OA 
7 a (x D) 
Therefore, 
wo, = [(xLowp) +(xDo,)']'” 


4 _ [(20\2, (@u\rhn 
A D L 

0.05 oo 22 
— = 90.061 or 6.1% 


W, = 9.43529E-—7m? 
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B. UNCERTAINTY IN POWER 


Wyr9 = 0.008% + 8uV [Ref 13, pg.587] 


Oy PwWre = 0.008% +8uV [Ref 13, pg.587] 


Wprq9 = +£0.001Q max value 
Q = Verwire X | Prwire 
V9 


= xX 
Y PrWire R 
22 











OQ. 29 
OV p, Wire Rio 
0Q _ Vp, Wire 
OV 0 Rio 
0Q Vo0 
= -.V +—- 
a Moe PrWire Roe 


V20 2 Verwi 2 V20 2 {1/2 
OO == ( Re x Ovni) or ( Sate X Wyo} + | -Voerwire x Rn x WI9 : 
Oa ( W y,P1Wire ) +( MY 22 r+ ine ) 1/2 
Q V prwire Vio Roo 


The worst case of uncertainty will occur when the voltages measured for Vp,y,, and Vo 
are small, less than 0.3 Volts. 


0. = 0 0/. \2 0 oy, \2 wnene:0O1 oy, 2 [12 
a (.008% +-0.0027%) + (0.008% + 0.0027%)" + ( =3564-700 


= 0.0524% or 5.2384E — 4 
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C. UNCERTAINTY IN HEAT FLUX 


> ¢ 
q A 
oo 
dQ A 
moe. 2sO 
OA 2 


[(0.061)° a (5.2384E-4)* |! 


6.1002E-2 or % 6.1 


D. UNCERTAINTY IN TEMPERATURE 


The uncertainty in the thermocouple measurements (w7-) 1s determined by the un- 
certainty in the platinum wire temperature probe (@>,p,,,,) and the error in the 
thermocouple’s in matching the probe temperature (@,7¢). 


te. = (GE @ ae 
= 0.4001°C 


E. UNCERTAINTY OF WIRE SURFACE TEMPERATURE 
From the calibration formula 


as 


oO, = 0.15 





Oro = 0.4001°C 
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1, = éERoprwire + Tp 


R 
= EV pr wire V+ Tot oT 














Where: 

€ = slope of calibration curve 
Tyo = axis intercept for R = 0 
Verwire = Voltage drop across Platinum Wire 
R39 = Resistance of in series Resistor 
V.q = Voltage drop across series Resistor 
6T = Difference between Data Acquisition and Final Calibration for a particular data po 
OT sur = V privire R 

de We 2 

OT suv ~~ 

hve & 
GN a 

0oT 

OT surf =e Roo 
OV piwire Voa 
OV surf = £ V piwire 
OR V2 
OT sung EY, Roo 
6NGn PrWire V2, 


1 


V prwireRra Roe 
OTsuf = Sy, * + (1 @79)° + ET Ov piwire “+(e 
Vie V0 


Rio 2 41/2 
(Vr a a oven) + (1 57) 
Vie 


Using data from Wire | at the following point, 
Woe. = “0.71673 Volts 
Vig = 0.24529 Volts 


Rxg = 1.9994 Ohms for wire 1 
bee = 993.902 C/O 
6T = 0.02°C 


This results in an uncertainty in temperature of, 


O Tsurfwire) = 0.854°C 


73 


Vp, Wire ) 
222 
Ve 


APPENDIX D. DATA ACQUISITION PROGRAMS USED 


A. PLATINUM WIRE RESISTANCE CALIBRATION PROGRAM 
The following program was used to obtain the Resistance value’s for the platinum 


wire specimans used in chapter 2. This program is written in HP basic language for use 
with controlling the HP-3852A Data Acquisition unit. 


1@ 
20 


3@ 
4Q 
50 
6Q@ 
70 
8 
92 
18Q 
11@ 
es 
13@ 
142 
152 
16@ 
17 
18@ 
19@ 
20@ 
21@ 
220 
iG: 
240 
25Q 
260 
2(Q@ 
28Q 
E30 
500 
51@ 
Oe 
350 
540 
55@ 
36@ 
370 
38@ 
B20 
4202 
410 


DIM Ohms(@:61 ) 
CUTE UT eGo NEG Erle 
OUTPUT 789;"RST @@0" 
OUTPUT 789; "USE 200" 
OUTPUT (2S. CONE ORG. 
OUTPUT 09s NERC wie 
OUTPUT 7033 NRbGsscee 
OUTPUT 7@3;"DEEAY Ga 
OUTPUT 709;"RANGE 12" 
OU TGs Wie Se = es 
OUTPUT 0 Ss" 2 = seal 
OUTRUN 7O9= “Mens sOnhibeles 
ENTER 789;:0hms(*) 
Alo=1.E+5@ 

= es 

Ahi=1.E-3@ 

A=@. 

FOR K=1 TO 6@ 
A=A+Ohms(k ) 

IF Ohms( kK )<Alo THEN 
Alo=Ohms(K ) 

END IF 

IF Ohms(K )>Ahi THEN 
Ahi=Ohms(K ) 

END IF 

NEXT K 

B=A/6@. 

PRINT "WIRE SAMPLE "3J 
RUN mmoe RESISTANCE AVERMGE ous 
PRINT EOP RES ScIPNGe Tien oURe Owl seme 
PRINT “HI RESISTANCE MEASURED IS "Ahi 
tPF J=5. scHeEN 

GOTO 362 

END Se 

GOTO 11@ 

GOTO 370 

DIM Temps(@:61) 

K=Q. 

OUTPUTS Oss oleC oul 
GUIPUT 709. USE 00g. 
OUTEU) 709s INTEGER M 
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420 
430 
44Q 
450 
460 
47@ 
480 
499 
50@ 
512 
520 
530 
540 
5SQ 
S60 
S72 
580 
590 
502 
612 
620 


OUIERE 2095 CONE TENET 
CUTeUT VOss NrEG=).0- 
OUTPUT 4093 NROGS 62” 
OUTPUT (09s WEEAY <O7< 
OUTPUL 4095: RANGE AUTO: 
OUTPUT. 704-4 —290 © 

QUT EU 7 Ogee 
OUTRO = (O35 MERsotlent hl. f 
ENTER 709; Temps( *) 
K=K+ | 

A=0. 

FOR L=1 TO 60 
A=A+Temps(L) 

NEXT 

Temp=A/6@. 

PRINT "THERMOCOUPLE";K; “TEMPERATURE”; Temp 
IF K=4 THEN 

GOTO 622 

ENE se 

GOTO 48Q 

END 


ie 


B. MAIN DATA ACQUISITION PROGRAM 


The following program was used to obtain the Temperature value’s for the platinum 


Wire specimans used in chapter 3. This program is written in HP basic language for use 


with controlling the HP-3852A Data Acquisition unit. 


12 

-@ 

5Q 

4@ 

5@ 

6@ 
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8@ 

9@ 

10 
11@ 
12@ 
15@ 
142 
15@ 
16@ 
17@ 
18@ 
190 
2Q@ 
21@ 
ae 

250 
240 
cod 
-60 
27Q 
28Q 
290 
5@@ 
319 
520 
550 
549 
550 
56@ 
570 
58@ 
390 
4Q@ 
410 
420 
450 
449 
459 


DIM Temos(@:16),Tk(1:7) ,7Tt(1:100) Flux (1:100) ,Tb¢(1@@) 


L=@ 

J=@ 

L=L+1 

IF L=81 THEN 

6010 1456 

Nii 

OUTPUT 7055 -FS1 eG * 
OUTPUT 7@9;"USE @@0" 
CUTEUT 70s... ENT EOE tae. 
GCUTPUT 7055 — CONE MI eENrT 
OUTPUT 7033 NPEC Ne. 
OUTPUT TOSS NREES. 17: 
OUNEUT (O39) DEUri =O. 
OUTPUT “(O03 RANGE AUG] 
OUNP UTS 70S. 21—244" 
OUTER 70s T= T+ 1 
OUTPUTS 70s; “MEAS @TEHe] =r 
ENTER 7@9;Temps( *) 
J=J+1 

A=@. 

FOR K=1 TO 15 
A=At+Temps( kK ) 

NEXT K 

Th( J =A/15. 

|S Saleye sia 

GOTO 300 

ENO LF 

GOTO 17@ 

DIM Volits(@: 11> VEC. 1a) 
J=0. 

Ct=0. 

OUTPUT 7@9;"RST O0Q" 
OUTPUT 7@39;"USE 90" 
GUTRUT 7635) INTEGER Mi. 
CUTIRUT 705; CONS Dey: 
OUTPUT 7@93;"NPLC 1.0° 
OUTPUT 7 OSs NRBGowt 2. 
OUTEUT 703; 4CDEEAY Vale 
CUTRUT 709s) RaNGe —iumem 
OUTPUTS (Ost =Jo 
OUTPUT) 70d. iT —t4+1 
GUIFUT 705s Memo OCY a. 


ENTER 7093;Volts(*) 


J=J+1 
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46 
470 
480 
49 
58@ 
51 
52 
55@ 
540 
550 
560 
579 
580 
592 
68@ 
612 
629 
65@ 
64@ 
652 
660 
670 
68@ 
69@ 
780 
710 
120 
750 
749 
750 
76@ 
77Q 
780 
73@ 
800 
810 
820 
830 
84@ 
859 
862 
878 
86@ 
890 
902 
312 
920 
9350 
340 
952 
36@ 
37@ 
980 


IF J=S5 THEN 
GOTO 420 

END IF 

C-Ctt 

A=. 

FOR K=1 TO 1@ 
A=A+Volts(K ) 

NEXT K 

ULC Ct I=A/ 10. 

IF J=6 THEN 

GOTO 590 

END IF 

GOTO 420 

DIM Voltw(@:11) 

t=O. 

ct=s. : 

OUTPUT 709;"RST 220" 
OUTPUT 709;"USE 200" 
OUTPUT 7@9;"INTEGER I" 
OUTPUT 7@9;"CONF DCV" 
OUTPUT 7@9;"NPLC 1.0" 
OUTPUT 7@9;"NRDGS 12” 
OUTPUT 7@9;"DELAY @.1" 
OUTPUT 709;"“RANGE AUTO" 
OUTPUT 7@9;"I=109" 
OUTPUT 709;"I=I+1" 
OUTPUT 7@9;"MEAS DCV,I" 
ENTER 709;Voltwi{*) 
J=J+1 

IF J=S THEN 

GOTO 712 

END IF 

Ct=Ctt i 

A=Q. 

FOR K=1 TO 10 
A=A+Voltw(K ) 

NEXT kK 

VEC Ct )=A/10. 

IF J=6 THEN 

GOTO 88@ 

END IF 

GOTO 710 
Bt=(TKC1 )tTkKC2 FTK(3 )4TK(4) 7/4. 
V1=Vk (6 ) 
T1=Vk(1)/1.9994 
V2=Vk(7) 
T2=Vk(2)/2.@Q76 
VZ=Vk(8) 
I3=Vk(3)/1.9997 


V4=VkR(9) 


T4=Vk(4)/2.9018 
VO=VR IO? 
T5=Vk(5)/2.8@25 
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390 ee de] 

1@@Q@ P2=I2*V2 

1@1Q@ PS3=15*VS 

1@2@ P4=14*V4 

105@ VEs=aiS+V5 

1040 Ri=V1/T1 

1050 R2=V2Z/12 

1060 !RS5=V3/15 

1070 '!R4=V4/14 

(Ge) SS RS=US7 rs 

1030 Fluxl=P ly i SoeSze35 [a 

1100 Flux2=sP27 lf) sSeace4e—s 

110) Pluxs=-sy1e 371 000SE—5 

1120 Flux4=P4/1.546Q565E-5 

1150 > Fluxs=Fs/ is /1000sEes 

1140 $1=€s5. 3021 255 Ril 25 4745955 

115Q@ T2=(55.666456+*R2 )-253.00959 

1160 T3=(55.141047#R3 )-252.7291351 
1170 14=(55.483646+#R4 )-252 .848797 
1180 [5=(55.085508*RS )-260.04883 


D130) Te=(TkheS ak Gea ce 


1200 PRINT “ TEMPER GIURE oe 

1210 PRINT “ AMBIENT AIR BUCK rEuiy CONDENSOR SURFACE" 

ae PRINT SckG 7) ote 

ipa PRU “sis 
1240 PRINT “WIRE _ TEMPERATURE CURRENT VOLT saa 

1250 eR ee (AMPS ) (VOLTS) ( DEGREE ae 
1260. “PRINT Seder | ea 


Per OSP eR TN eure ean yie 

Lacd: OUP R Nie ee ooo 
129@ '!PRINT “4" 74,14 ,V4 
1590" PRIN Ss Sloe vS 


[S70> PRINT = BEAT Pex POWER RESISTANGHas 
1320 '!PRINT “ (WATTS) (OHMS ) (WAT Is7 iia 
1330) PRINTS tae user een 


1340-7 PRINT S22 sr tox 2 eo oRe 
155@. |PREINT 3" pee eS. 85 
1360 '!PRINT “4° .Flux4 P4 R4 
1570 (!PRINT °S" .-lux5 PS. RS 
IS8Q° (PRIN ie ee 
Is990" Tithoah 
1400 Tb(L)=Bt 
1410 Flux(Ll)=Fluxi 
1420 WAIT 102 
1430 GOTO 3@ 
1440 GOTO 1530 
1450 M=1 
1460 PRINT "“ Wi RES Teme HEAT FLUX BULK TEMP" 
1470 IF M=3 THEN 
1480 GOTO 153@ 
1490 END IF 
1590 PRINT OM TiC) eax CMe Span 
1510 M=M+1 
1520 GOTO 1470 
1530 END 
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C. NUSSELT EVALUATION PROGRAM 
The following program was used to obtain the Nusselt Number for the platinum 
wire specimans used in chapter 3. This program is written in Fortran language for use 


with controlling the IBM-3803 Main Frame Computer. 


FILE: NUSSELT FORTRAN Al 


SUBROUTINE NUSSELT (€X1,T1,M,Q) 
THIS SUBROUTINE WILL CALCULATE THE HEAT TRANSFER COEFFICIENT FOR THE 
DATA OBTAINED FROM THE EXPERIMENTAL RUNS BY USING THE CHURCHILL AND 
CHU CORRELATION 
DIMENSION X1(%),T10%),Q(%) 
REALX4 X1,71,Q,K,NU,NUD 
OPENC29,FILE="RALY*,STATUS=*UNKNOWN® , FORM=* FORMATTED ® 
DIAMETER OF THE WIRE BEING ANALYSED IS NOMINALLY .05 MM 
D=>.0E=5 
GRAVITY TERM OF 9.81 M/S2 
G=9.81 
DO 100 I=1,™ 
FILM TEMPERATURE OF THE PHYSICAL PROPERTIES 


CC TF CXECTO ett 1) ot 56.20)72. 
Tr CALCIO te. eC 1 72: 


TEMPERATURE DIFFERENCE BETWEEN THE WIRE AND BULK FLUID 


ce POS XTC Ss Vici) + 56.03 
TD = X11) 


THERMAL CONDUCTIVITY OF FLUID W/MK 

K = (0.6055 - 0.00115XTF)/10. 
BULK MODULUS OF THE FLUID 1/K 

BETA = 0.00261/(1.740 - 0.00261XTF) 
KINEMATIC VISCOSITY OF THE FLUID M2/s 

NU = 1.203952E-8 * EXP(1058.4109/(TF+273.15)) 
SPECIFIC HEAT OF FLUID 

CP = (0.24111 + 3.70374E-4 X TF)¥4186.00 
FLUID DENSITY KG/M3 

RHO = (1.740 - 0.00261*TF}*1000.0 


NN00N 


MO O730O0 9090 A090 A090 ANN ANAND ANNA A ANN ANN 


vB, 


ann MONON ANN AN 


ONaN 


THEMAL DIFFUSIVITY 
ALFA = K/CRHOXCP) 

PRANDLT NUMBER OF FLUID 
PR = NU/ALFA 

RALYEIGH NUMBER OF FLUID MUST BE 10%X-5< RA < 10%*12 


RA = (CG¥BETAXTDXD¥X3.)/CNUXALFA) 
PRINT X®, RA 


NUSSELT NUMBER OF FLUID NU=HXD/K 
TERM = €€0.518*CRAXX(0.25))¥€14+00.559/PR) XO .6)*¥(-5./12. ) ) X15 
1+CO. 1 ¥RAX¥C1.73.))¥*¥15.3**C1.715.) 


nea es 2./TERM) ) 


THERORICAL HEAT FLUX WATTS/M2 


QCI) = KXNUDXTD/D 


Pas WRITEC29,500) TD,RA,NUD,PR, Q(I) 


5900 FORMAT (5(F12.4,2x)) 


CLOSE (29) 
RETURN 
END 
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